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ABSTRACT 

Several concepts now exist for small, space-based missions to directly char¬ 
acterize exoplanets in reflected light. While studies have been performed that 
investigate the potential detection yields of such missions, little work has been 
done to understand how instrumental and astrophysical parameters will affect 
the ability of these missions to obtain spectra that are useful for characterizing 
their planetary targets. Here, we develop an instrument noise model suitable 
for studying the spectral characterization potential of a coronagraph-equipped, 
space-based telescope. We adopt a baseline set of telescope and instrument pa¬ 
rameters appropriate for near-future planned missions like WFIRST-AFTA, in¬ 
cluding a 2 m diameter primary aperture, an operational wavelength range of 
0.4-1.0 /rm, and an instrument spectral resolution of A/AA = 70, and apply our 
baseline model to a variety of spectral models of different planet types, including 
Earth twins, Jupiter twins, and warm and cool Jupiters and Neptunes. With 
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our exoplanet spectral models, we explore wavelength-dependent planet-star flux 
ratios for main sequence stars of various effective temperatures, and discuss how 
coronagraph inner and outer working angle constraints will influence the potential 
to study different types of planets. For planets most favorable to spectroscopic 
characterization—cool Jupiters and Neptunes as well as nearby super-Earths— 
we study the integration times required to achieve moderate signal-to-noise ratio 
spectra. We also explore the sensitivity of the integration times required to ei¬ 
ther detect the bottom or presence of key absorption bands (for methane, water 
vapor, and molecular oxygen) to coronagraph raw contrast performance, exozo- 
diacal light levels, and the distance to the planetary system. Decreasing detector 
quantum efficiency at longer visible wavelengths makes the detection of water 
vapor in the atmospheres of Earth-like planets extremely challenging, and also 
hinders detections of the 0.89 pm methane band. Additionally, most modeled 
observations have noise dominated by dark current, indicating that improving 
CCD performance could substantially drive down requisite integration times. Fi¬ 
nally, we briefly discuss the extension of our models to a more distant future 
Large UV-Optical-InfraRed (LUVOIR) mission. 


1. Introduction 


Exoplanet atmospheric characterization is a rapidly progressing field, and will likely 
continue on this trajectory as the detection and study of nearby exoplanets was highlighted 
as one of three main science objectives in the 2010 National Research Council decadal survey 
of astronomy and astrophysics^. Following the first detection of an exoplanet atmosphere 


(Charbonneau et ah 2002), secondary eclipse and transit observations have been used to 


characterize the atmospheres of a number of hot Jupiters (Grillmair et ah 2008, Swain et ah 


2008, Pont et ah 2008, Swain et ah 2009, Sing et ah 2009, Madhusudhan & Seager 2009). 


More recently, observations have begun to probe new categories of planets, including mini- 


Neptune and super-Earth exoplanets (Stevenson et ah 2010; Bean et ah 2010, Kreidberg 


et ah 2014; Knutson et ah 2014a||b Ehrenreich et al. 2014, Fraine et ah 2014). Of course, 
and as with any young field, some findings remain controversial or have undergone substantial 


revision (Line et al. 2014, Hansen et al. 2014, Diamond-Lowe et al. 2014) 


Anticipating future exoplanet direct detection missions, several studies have investigated 


the detection capabilities of coronagraph or starshade instruments. Agol (2007) presented 
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analytic models of the detection efficiency of coronagraph searches and investigated depen¬ 
dencies on key parameters, including survey duration, telescope diameter, and telescope 


inner working angle (IWA). Brown & Soummer (2010) discussed methods for optimizing 


bull et ah 

2012 

) have emphasized techniques for evaluating the effectiveness of proposed 

missions. ’ 

Vlawet et al. 

(2014 

) studied the statistics of high-contrast imaging at small angles. 


and found that the limited number of resolution elements near to the star (i.e., at < 2X/D) 
can lead to large errors in contrast estimations and false alarm probabilities computed by less 


sophisticated methods. Stark et ah (2014) presented an improved technique for optimizing 


exoplanet detection completeness, and explored the dependencies of exo-Earth yield on a 


number of astrophysical, mission, and instrument parameters. In a follow-up paper. Stark 


et ah (2015) used their completeness model to compute lower limits on the aperture size 


required to hnd and characterize a signihcant number of Earth-twins. Very recently, Greco 


& Burrows (2015) studied the detection capabilities of a coronagraph-equipped WFIRST- 


AFTA mission^ (Spergel et ah 2013). 


Technologies and ideas for direct detection and characterization of exoplanets have long 
been discussed, and there exists a body of white papers and summary documents that 


outline strategies for direct characterization missions (Beichman et ah 1999, Des Marais 


et al.||2002 i Schneider et al.||2008 i Lunine et al.||2008 Cockell et al.||2009 ). Recently, following 
a growing interest in exoplanet direct detection^’^, NASA’s Exoplanet Exploration Program 
commissioned a number of “quick study” reports focused on exoplanet characterization at 


visible wavelengths. Marley et ah (2014) discussed how reflected light observations could 


be used to characterize existing radial velocity-detected exoplanets. Hu (2014) studied how 


direct imaging could constrain methane abundances and cloud properties using low resolution 


spectra. Finally, Burrows (2014) explored the diversity of giant exoplanet spectra that may 


be accessible to a coronagraph-equiped WFIRST-AFTA mission. 

Only a subset of the aforementioned quick studies considered the influence of noise on 
spectral observations, which, if included, is commonly applied by assuming a constant signal- 
to-noise ratio (SNR) across the entire modeled “observation”. While some studies have done 
extensive modeling of the influence of speckle noise (and techniques for minimizing such 


noise 


on direct exoplanet observations (Marois et ah 2000, Sparks & Ford 2002), few pub¬ 


lished studies exist that comprehensively examine how spectral observations are influenced by 
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stellar leakage, zodiacal light, exozodiacal light, and other noise sonrces. Maire et al. (2012) 
investigated how noise would influence the characterization capabilities of the proposed the 
Spectro-Polarinietric Imaging and Characterization of Exoplanetary Systems (SPICES) mis¬ 
sion (Boccaletti et al. 2012). This work included noise due to speckles, exozodiacal light, 
zodiacal light, and read-out, resulting in a handful of spectral simulations for planet-star 
systems for a specified distance and integration time. Recently, Leger et al. (2015) studied 
design missions for a 2-meter class space-based coronagraph mission and a four-spacecraft 
formation flying interferometer mission, assuming noise from only stellar leakage and exo¬ 
zodiacal light, with a primary goal of determining how effective such missions would be in 
studying exo-Earths given new estimates of 


Here, we build on the limited previous studies of the influence of noise on the ability of 
space-based coronagraphs to characterize exoplanets. We begin by outlining a framework for 
modeling noise for such missions, including stellar leakage, exozodiacal light, zodiacal light, 
dark current, and read noise. We then investigate the capabilities of a 2-meter class mission 
to observe and study a variety of Jupiter- and Neptune-sized worlds as well as terrestrial exo¬ 
planets. Using a characteristic set of mission parameters, we show how expected integration 
times will scale with wavelength and host star spectral type, and, for key spectral features, 
we show the sensitivity of these integration times to distance to the star-planet system, exo¬ 
zodiacal light levels, and coronagraph raw contrast performance. Finally, we discuss how our 
framework extends to the characterization of Earth twins with Large UV-Optical-InfraRed 
(LUVOIR) telescope missions. 


2. Model Description 


We consider a space-based telescope of diameter D designed to detect and characterize 
exoplanets in reflected light. The telescope is equipped with a coronagraph that achieves a 
raw contrast C (which can, in general, be wavelength-dependent) and has inner and outer 
working angles of 6*iwa and 6*owa, respectively. The planetary light is is distributed into 
a point spread function (PSF), where an aperture of width (or diameter) XX/D is used 
in the image plane to extract the planetary signal. However, a number of other sources 
will contribute to the detector counts within the defined aperture, including Solar System 
zodiacal light, exozodiacal light, leaked stellar light, dark current, and read noise. 


The Appendix contains a presentation of our models of the planet count rate and count 
rates from key noise sources, which are similar to those of Brown (2005). For convenience, a 
synopsis of all of the model parameters, variables, and outputs are given in Table 1. Briefly, 
our planet count rate (cp), zodiacal light count rate (c^), exozodiacal light count rate (cez). 
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stellar leakage count rate (cik), dark current count rate (cd), and read noise count rate (cr) 
are computed using Equations |A5[ |A7[ |A10[ |A11[ |A12[ and A21 respectively. 


2.1. Baseline Parameters and Input Data 


It is useful to define a baseline set of astrophysical, telescope, and instrument parameters, 
though later sections of this paper will explore sensitivity to key variables (e.g., distance, 
exozodiacal light levels, and coronagraph raw contrast). Baseline astrophysical parameters 
are given in Table 2, while telescope and instrument parameters are in Table 3. 


Regarding astrophysical parameters, we adopt the standard, empirically-derived solar 
spectrum of Wehrli (1985). As mentioned in Section A.2, we use a surface brightness of 
Mz V = 23 mag arcsec“^ for Solar System zodiacal light. For exozodiacal light, we assume 
one zodi (i.e., = 1) with a surface brightness of = 22 mag arcsec“^, which is a 

factor of ~2 larger than Solar System zodiacal surface brightness, as exozodiacal dust both 
above and below the midplane contribute (Stark et ah 2014). 


For baseline telescope and instrument parameters, we consider a 2-meter class telescope 
equipped with a coronagraph designed to achieve a raw contrast of C = 10“® between 
6*iwa = 2A/Z1 and 6*owa = lOA/D. The total system throughput is T = 0.05, and we 
mimic the decreasing sensitivity of CCDs at red wavelengths by parameterizing the quantum 
efficiency according to 


I u.a, A ^ u. / 

^“\o.9(l-^), 0.7/rm 


( 1 ) 


< A 




We use a characteristic dark current of 5 x 10“^ s“^ and read noise of 0.1. We assume a square 
aperture width of 1.5A/D (i.e., X = 1.5), which would be a sensible choice to help eliminate 
sensitivity to the sub-pixel (or sub-lenslet) location of the planet, whose PSF is sampled at 
0.5A/D. Finally, we assume that the coronagraph is paired with a spectrometer that covers 
0.4-1.0 /xm at 77 = 70, which is sufficient to resolve many molecular absorption bands and 
features throughout the visible wavelength range (see Figure]^. Note that, while a resolution 
of 70 is small enough to resolve the A-band of molecular oxygen, this may not be the ideal 
resolution for detecting this band ( Brandt &: Spiegel|2014 ). Overall, our assumed parameters 
are meant to be representative of either the previously-mentioned WFIRST-AFTA mission 
or the Exo-C mission concept (see reports linked in the Introduction). 

Our reflectivity data come from a number of different sources. For Earth twins we use the 
extensively validated Virtual Planetary Laboratory 3-D spectral Earth model from Robinson 
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et al. (2011), which realistically simulates wavelength- and phase-dependent spectra of the 
Pale Blue Dot (Robinson et al. 2010, 2014). In certain cases, we investigate super-Earths, 


which we take to have an Earth-like reflectivity and a radius of I.SRq, and Venus twins, with 


a characteristic reflectivity from the validated models of Arney et al. (2014). For Jupiter 


and Neptune twins, we use the observed geometric albedos from Karkoschka (1998) and a 
Lambertian phase function given by 


‘hida') = 


sin a + {n — a) cos a 


TT 


( 2 ) 


although realistic planetary phase functions are not perfectly Lambertian (e.g., Cahoy et al. 


2010, their Figure 15). Recall that the phase function captures phase-dependent brightness 


changes as the planet moves through its orbit and that the brightness increase from quadra¬ 
ture to full phase is a factor of ~3. Finally, we also consider “warm” and “cooF’ Jupiters 


and Neptunes, using the wavelength- and phase-dependent models of Cahoy et ah (2010) for 


both a 0.8 AU and 2 AU (from a Sun-like star) Jupiter and Neptune (at 3x and 30 x heavy 
element enhancement, respectively). All planets are taken to be at a characteristic phase 
angle of a = 7r/2 = 90° (i.e., quadrature). 


2.2. Noise 


The total number of counts, Ctot, is computed from the total count rate using 

Ctot (Cp T Cz -)- Cez T Qk T C£) J- Cj^) Afexp (Cp “F Cb) Afexp ; (3) 

where Atexp is the exposure time, and we have defined the total background count rate, 
Cb = Cz + Cez + cik + cd + cr. The total number of background counts would be Cb = CbAt 


Following Brown (2005), the noise counts are taken as 

Cnoise = y/Cp J- 2C'b , 


exp‘ 


( 4 ) 


which assumes a background subtraction (e.g., via a spacecraft roll maneuver) is performed 
to detect the planet (implying that the single exposure photon counting noise limit is not 
achieved), and Cp is the total planet counts. With the signal-to-noise ratio given by 


SNR= 


C ■ 

noise 

we can compute the exposure time needed to achieve a given SNR by 

Atexp = ^i^^^SNR^ . 


( 5 ) 


( 6 ) 
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Note that the exposure time required to achieve a given SNR is distinct from the SNR 
required to detect a certain molecule (e.g., Brandt fc Spiegel]2014 ), or to distinguish between 
models with different heavy element enhancements (Maire et ah 2012). The former (given 
by Equation]^, when applied at the bottom of several molecular features (as is done in later 
sections), is representative of the requisite exposure time for determining molecular abun¬ 
dances and cloud properties for exoplanet atmosphere (Marley et al. 2014). Such determi¬ 
nations require comparisons to cloudy atmospheric models, both at continuum wavelengths 
and across/in molecular bands, thereby simultaneously retrieving cloud and atmospheric 
parameters (including the abundances of key species). 


We also explore the integration time required to simply detect a molecular feature by 
the deviation it causes from a flat continuum, which is sufficient for stating whether or not a 
certain species is present in an atmosphere. Following Misra et al. (2014), we write the band 
SNR by dehning the signal as the count difference between the spectral model and a model 
without absorption, determined by htting the continuum on both sides of the absorption 
band. The noise counts are summed over the band. Thus, if Ccont is the continuum count 
rate, then the band SNR is given by 


SNRhanri — 


Ccont 

3 


''PJ 


Cpj -|- 2Cb, 

3 


-^''exp 


(7) 


where the sum is over all spectral elements (denoted by sub-script ‘j’) within the molecular 
band. Band detection times, in general, will be long for shallow and narrow features, whereas 
the previously discussed exposure times to reach a given SNR at the bottom of a band is 
shortest for shallow features. These two times will be roughly equal for a feature that is one 
resolution element in width and that is moderately deep. 


3. Results 

The following subsections explore the detectability of, and integration time required to 
detect, key spectral features for various categories of planets. We dehne a planet “type” as 
having hxed radius, reflectivity [A$(a)], and incident flux from its host star. The latter 
requirement is sometimes referred to as the “flux (or insolation) equivalent distance”. Note 
that we focus on the visible wavelength range (0.4-1.0 /rm), as exoplanet characterization in 
reflected light at near-infrared wavelengths is especially difficult for small telescopes due to 
IWA constraints. 

To get a sense for the scales involved, we note that an Earth twin around a Sun-like 














star, using A = 0.2, <l>(7r/2) ~ I/tt, Rp = R^, and r = 1 AU, would have the canonical 
planet-star flux ratio (from Equation |A3[) of Then, assuming our baseline parameters. 


5 X 10"^ s"^ 


and using A = 0.55 /um with a 8 nm-wide bandpass. Equation A5 gives Cp 
(or about 2 hr“^) for an Earth twin at 5 pc. By comparison, the leakage count rate is 
20 hr“^ for C = 10“®, while the zodiacal and exozodiacal light count rates are 2 hr“^ and 
4 hr“^, respectively. With a dark count rate of about 240 hr“^ and a read noise count rate 
of 10 hr~^, the exposure time needed to reach SNR = 5 is roughly 5 x 10^ hr—an unfeasibly 
long integration time, but which could be reduced by decreasing spectral resolution. For a 
Jupiter twin at 5 pc, with a geometric albedo of 0.5, the planet-star flux ratio increases to 
~10“®, and the planet count rate increases to 20 hr“^ (although this twin, at 5.2 AU from its 
Sun-like host, would be outside 6*owa = lOA/iA at quadrature) while the exozodiacal light 
count rate decreases to 0.2 hr“^. Here, then, the time needed to reach SNR = 5 is only 30 
hours. 

As demonstrated by these count rate example calculations, the dominant term in the 
noise counts (Equation —for a wide range of planet types, distances, stellar effective tem¬ 
peratures, and wavelengths—is dark current. This term is strongly influenced by two pa¬ 
rameters: the dark current rate {D^-) and the width of our assumed photometric aperture 
(controlled by X). The latter term determines the area used to extract the planetary signal, 
which leads to a dependence in the dark noise count rate. For our assumed baseline 
set of parameters and a characteristic distance of 10 pc, all of the planet types investigated 


in Section 3.2 and later have dark current as the dominant noise source, leading the expo¬ 


sure times required to achieve a given SNR to scale with the planet count rate as (roughly) 
cn/Cp. Note that this scaling leads to a strong dependence on planetary properties, implying 
that, for example, observing planets closer to full phase would drive down integration times 
by as much as an order of magnitude from the estimates provided below (which are given 
at quadrature phase). Also, the overall importance of dark current in our results below 
is influenced by our baseline throughput assumption (selected to be consistent with, e.g., 
WFIRST-AFTA)—improved telescope and instrument throughputs will drive up the plan¬ 
etary signal, while also increasing the relative weights of noise due to stellar leakage and 
zodiacal and exozodiacal light. 


3.1. Planet-Star Flux Ratios 

Insight into the complexities associated with directly observing cool exoplanets comes 
from investigating planet-star flux ratios for our different types of worlds. The first panel 
of Figure shows the planet-star flux ratio as color contours for a Jupiter-sized world with 
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A® (a) = O.S/tt (independent of wavelength) and with a Jupiter-like top-of-atmosphere 
flux (fixed at -Ftoa = 50.5 W ni“^), as a function of main sequence host star effective 
temperature. As we have assumed a gray reflectivity, the planet-star flux ratio only depends 
on Teff (assuming a simple relationship between this and Rg on the main sequence), with 

Fp _ Ad>(a)FTOA [RpV 
ctT^s \Rs) 

where the dependence on planetary orbital distance is functionally removed. Here, since 
the planet flux is fixed by (a) and FpoA (regardless of stellar effective temperature), the 
planet-star flux ratio increases for cooler stars as these stars are intrinsically fainter. 


The planet-star flux ratio contours for the gray case in Figure tell an incomplete story, 
however. First, while the results indicate that detecting a Jupiter-sized planet around cool 
stars (with < 4000 K) would be straightforward with a coronagraph capable of achieving 
a raw contrast of only 10“®, it is important to remember that the probability that such cool 


stars host a gas giant is small (of order a few percent, Johnson et al.| 2010 Montet et ah 


2014). Additionally, the coronagraph inner and outer working angles will constrain the range 


of planet-star apparent separations (which are related to orbital separations) accessible to 
the telescope when observing a system at a given distance. These constraints are shown as 
dashed lines for distances of 5 pc and 15 pc in Figure]^ For fixed FpoA and a system distance 
of 5 pc, the OWA would not allow observations of the gray Jupiter world at quadrature over 
the entire wavelength range for host stars with effective temperatures above 5800 K. When 
placed at 15 pc, the IWA can prevent observations at the longest wavelengths for stars with 
effective temperature below 4900 K. 


The remaining panels in Figure]^ show non-gray cases, where realistic reflectivities are 
used. Like the gray case, the Jupiter twins have Ftoa = 50.5 W m“^, while the cool and 
warm Jupiters have FpoA = 342 W m“^ and FpoA = 2140 W m“^, respectively. With the 
planet-star flux ratio now a function of wavelength, these plots allow us to easily determine 
roughly what raw contrast must be achieved to detect different types of Jupiters at continuum 
wavelengths and within absorption features. The Jupiter twins and the cool Jupiters are 
relatively reflective and have several methane bands covering a range of feature depths. 
Note that the Jupiter twins have low planet-star contrast ratios due to their small FpoA- 
Additionally, despite their large Rtoa, the warm Jupiters have low planet-star contrast ratios 
at longer wavelengths due to strong methane and water vapor bands as well as a lack of water 
clouds. 


Figure is similar to Figure but for Neptune twins, cool Neptunes, and warm Nep- 
tunes. Neptune twins have very low planet-star flux ratios due to a very small FpoA 
(1.5 W m“^), and would also be difficult to detect for a wide range of stellar effective 
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temperatures due to the telescope OWA restrictions. The flux ratios improve for the cool 
and warm cases, but are lower than the equivalent Jupiter-sized cases as Neptune is 2.9 times 
smaller than Jupiter. Also, the methane and water bands for these cases are deeper than 
those for the Jupiters, yielding flux ratios near over much of parameter space. 

Continuing downwards in the size/mass regime. Figure shows planet-star contrast 
ratios for Venus twins (iVoA = 2610 W m“^), as well as for Earth twins and super-Earths 
(both with Ftoa = 1370 W m“^). The Venus twins have comparatively high planet-star flux 
ratios, both due to their high reflectivity and large Ftoa, and have nearly featureless spectra 
over the highlighted wavelength range. Inner working angle constraints make Venus twins 
very difficult to detect for our chosen architecture, especially at longer wavelengths. For 
Earths and super-Earths, the visible wavelength range contains key spectral features—the 
0.76 /im O 2 A-band and the 0.94 /xm water band—but IWA constraints largely limit our 
ability to characterize such planets around cool stars. 

The relationships between wavelength-dependent reflectivity, host star effective temper¬ 
ature, and planet-star flux ratio discussed above can be straightforwardly manipulated to 
show contours of stellar effective temperature where planet types achieve a given flux ratio. 
Four such plots are shown in Figure which illustrate the range of different stellar spectral 
types that can be investigated for a specified planet type and contrast performance. These 
diagrams demonstrate how, for a fixed planet-star flux ratio, the depths of absorption bands 
can best be probed for cooler stars, where the low intrinsic stellar brightness compensates 
for the low planetary reflectivity. As before, inner and outer working angle limitations are 
shown for planet-star systems at 5 pc and 15 pc. 


3.2. Integration Times: Full Spectral Range 


We now shift our attention to the integration times required to achieve a specified SNR 
for different planet types. These exposure times determine whether or not it is realistically 
feasible to acquire the types of observations needed to characterize planets at a wide range 
of wavelengths, especially given constraints on mission duration and overall mission time 
allotted for exoplanet studies. While we highlight SNR = 5 as a minimum value needed for 
characterization—Marley et ah (2014) were able to draw scientifically valuable conclusions 
from their simulated data at this SNR—the exposure times discussed here can be easily 
scaled to other SNRs since Afexp oc SNR^. 

Wavelength-dependent integration times required to achieve SNR = 5 fixed iVoA cases 
are shown for Jupiter twins and cool Jupiters at a characteristic distance of 10 pc in Figure 
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Requisite integration times for the cool Jupiters are promising, reaching SNR = 5 in 10- 
100 hr, even in the depths of some methane absorption bands. For the Jupiter twins, 
detecting the bottom of these bands requires more than 10^ hours of observation, owing 
to these planets’ larger orbital distances, although binning data to lower spectral resolution 
[TZ =20-30) would drive down the requisite exposure times (by a factor of a few). 

Figure shows integration times required to achieve SNR = 5 for cool Neptunes at a 
distance of 10 pc. Note that these integration times show a much larger dynamic range than 
the cool Jupiter cases owing to the generally deeper methane bands in the geometric albedo 
spectrum of the Neptunes . Overall, integration times shorter than 10^ hr are only achieved 
below 0.7 fim, owing to both strong absorption features and falling stellar spectra at longer 
wavelengths. 

Finally, Figure]^ shows integration times to SNR = 5 for Earth twins and super-Earths, 
both taken to be at 5 pc distance. To maintain integration times shorter than about 10^ hr, 
the Earth and super-Earth observations have been degraded to spectral resolutions of 10 
and 20, respectively (which would not allow detection of the 0.76 /im A-band of molecular 
oxygen). Outside of the 0.94 /im water band, the Earth twins have integration times of order 
10^ hr, although achieving such observations for Sun-like stars and hotter would require a 
contrast noise floor better than 10“^°. By comparison, super-Earths avoid the contrast 
limitation over most of parameter space, and require integration times less than 500 hr for 
most wavelengths (at 77 = 20). 


3.3. Integration Times: Molecular Features 


By highlighting key molecular bands in certain planet type spectra, we can investigate 
the integration times required to detect these features. We explore both the integration 
time required to achieve a given SNR at the bottom of a molecular band (Equation]^, and 
the integration time required to detect a band (or bands) of certain species (Equation]^. 
Recall that the former is indicative of the integration time required to achieve observations 
suitable for abundance determinations (where detections of the bottoms of multiple bands 
are required), while the latter indicates the time required to determine whether or not a 
species is present. As before, we set the detection SNR to be 5, which can easily be scaled to 
other SNR values. Also, for all the bands discussed below, our baseline spectral resolution 


(77 = 70) is sufficient to resolve these bands (Des Marais et ah 2002). 


Figure]^ shows integration times for detection of species and bottoms of bands for cool 
Jupiters as a function of distance to the planetary system. We focus on methane, which 
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is the dominant absorber for the cool Jupiters (and Neptunes) at visible wavelengths. For 
detecting the bottoms of bands, we highlight the 0.73 /xm and 0.89 /rm methane bands, which 
are, respectively, moderate and strong bands. For detecting methane, we use the 0.73 /xm 
band, although required integration times for detecting the 0.89 /xm band are similar, owing 
to a trade-off between detector quantum efficiency (lower at these wavelengths) and band 
strength (higher). We also show that the integration time required to detect methane is 
decreased by about 50% if both bands are used. Including the 0.62 /xm methane band does 
not further decrease integration times, as the shallowness of this band leads to long required 
integration times for detection. 

Figure is similar to Figure except for cool Neptunes. For detecting band bottoms, 
we use the 0.54 /xm and 0.62 /xm methane bands, as detecting the bottom of longer wave¬ 
length bands would require contrast ratios better than 10“^°. We show integration times for 
methane detection using the 0.62 /xm band, and including both this band and the 0.73 /xm 
decreases requisite times by about 50%. 

Finally, Figures [TT] and \U2 \ show integration times for band bottom and species detections 
for Earth twins and super-Earths, respectively. Here, we emphasize the 0.76 /xm A-band 
of molecular oxygen and the 0.94 /xm water vapor band. The A-band is a key potential 
biosignature, while the water vapor band is essential for discerning planetary habitability. 
Note that the time required to detect the bottom of the A-band is much shorter than the time 
to detect the presence of the band. This is due to the convolution of the high-resolution Earth 
reflectance spectrum with the spectrometer line-shape function, which serves to decrease the 
depth of this band at 77 = 70 and, thus, makes detection of the band difficult. Thus, the 
integration time required to reach the bottom of the A-band at a SNR of 5 is not much 
different from the time to reach the same SNR in the surrounding continuum. We note 
that requisite integration times for Earths or super-Earths placed at the outer edge of the 
habitable zone are typically 5-10 times longer due to the lower stellar flux incident on these 
planets. 


3.4. Sensitivity to Key Parameters 

Our baseline parameter set assumed that every planet-star system has one exozodi and 
that the coronagraph could deliver a raw contrast of 10“®. However, exozodi levels can 
vary by orders of magnitude, as can the designed raw contrast (depending on which mission 
architecture is adopted). Thus, we explore sensitivity to these two parameters by examining 
their influence on our ability to characterize a planet for key atmospheric constituents. 
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Figure shows the integration time required to achieve SNR = 5 in the bottom of 
the 0.73 /im methane band as well as the integration time required to detect this band 
SNRband = 5 for our cool Jupiter model. Panels demonstrate sensitivity to exozodi levels and 
raw contrast performance at a hxed distance of 10 pc. For these studies, the coronagraph IWA 
prevents observations of this planet type around stars cooler than 5,300 K. As the 0.73 /xm 
methane band is not particularly deep or wide in the cool Jupiter reflectance spectrum, 
integration times for band detection are longer than for detecting the bottom of the band. 
Nevertheless, even at poor raw contrast performance and/or larger {Nez ^ 10) levels of 
exozodi, integration times for either band bottom or species detection are below 100 hr. 


Figure 14 is similar to Figure 13, but for the cool Neptune model. Here we use the 
0.62 /im methane band, because, as mentioned earlier, the bottom of the 0.73 /xm band is at 
a contrast ratio of less than 10“^°. We place the planet-star systems at a distance of 5 pc. 
Integration times for band detection are now shorter than for detecting the bottom of the 
band, as the 0.62 /xm methane band is deeper in the cool Neptune reflectance spectrum than 
the 0.73 /xm band in the cool Jupiter spectrum. 


Finally, Figure shows sensitivity tests for the integration times required to detect 
the 0.76 /xm oxygen A-band for 1.5i?0 super-Earths. We place the planet-star systems at a 
distance of 3.7 pc, which is roughly the distance to r Ceti. We do not consider the 0.94 /xm 
water vapor band as Figure \U2 \ shows that detecting this band at 3.7 pc would require nearly 
a 10^ hr integration time. Also, we do not consider the integration time required to detect 
the bottom of the band as this is shorter than the time to detect the band. 


3.5. Towards Large UV/Optical/Infrared Telescopes 


Large UV-Optical-Infrared telescopes have been proposed as important tools for con¬ 
tinuing NASA’s vision for exploring and characterizing exoplanets beyond the James Webb 


Space Telescope mission (Kouveliotou et al. 2014). The capabilities of such a 10-meter class 


space-based telescope, if equipped with a coronagraph, can be investigated using the formal¬ 


ism outlined above. As an example. Figure 1^ shows integration times to SNR = 5 across 
a broad wavelength range, and sensitivity to exozodi levels for the 0.76 /xm oxygen A-band 
for Earth twins, as would be the primary target of a LUVOIR mission. We have adopted 
an extended wavelength range (0.4-3 /xm), and have assumed the distance to these systems 
is 10 pc. We have extended the OWA to a more generous 20A/iJ and use C = 10“^°, as 
could hopefully be achieved by such a distant-future mission. Also, we extend the IWA to 
3X/D, as telescope stability for a 10-meter class observatory would make achieving smaller 
IWAs very challenging. (Note that this larger IWA is outside the 2X/D “practical limit to 
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small IWA coronagraphy” determined by Mawet et al.[ (2014).) Finally, we assume that the 
visible and near-infrared wavelength ranges would be handled by two separate detectors, 
with respective pixel/lenslet sizes determined at 0.4 /xm and 1.0 /xm. For the near-infrared 
detector, we adopt a dark current, read noise, and quantum efficiency that are representative 
of HgCdTe detectors (Morgan & Siegler j2M5 ). Specifically, we use D^- = 10“^ s“^, R^- = 1, 
and q = 0.85. The visible wavelength detector is the same as our baseline model, except we 
adopt a more optimistic dark current of D^- = 10" 


s-k 


For the adopted LUVOIR architecture and for Earth twins at 10 pc, most of the 
wavelength-Teff parameter space is dominated by noise from dark current, which is not 
influenced by telescope diameter. If not for the larger dark current rates of the near- 
infrared HgCdTe detector, exozodiacal light would be the dominant noise source at longer 
wavelengths, owing to its strong wavelength dependence (oc A^Fs^a)- Note also that the 
wavelength-dependent IWA at 3X/D limits our ability to obtain a full spectrum through the 
near-infrared for all but the hottest stars. For Earth twins at 5 pc, the IWA restrictions 
are less strict (although the OWA can restrict observations at quadrature for stars hotter 
than 5,500 K), and leakage from the star can dominate over other noise sources at visible 
wavelengths for stars hotter than about 6,000 K. This all assumes that the telescope system 
can be cooled to a sufficiently low temperature to minimize thermal noise contributions to 
observations in the near-infrared (i.e., below about 80 K). As an example of the capabilities 


of a LUVOIR mission. Figure 17 shows a simulated observation of an Earth twin orbiting 


a solar twin, assuming the system is at 10 pc and an integration time of 200 hr. Note that 
IWA constraints cut off the observations at the longest wavelengths. 


Discussion 


Direct characterization of exoplanets in reflected light will require moving beyond V- 
band detections, and will utilize a wide range of wavelengths. This added spectral dimension 
influences our understanding of spectrographic observations with coronagraph-equipped tele¬ 
scopes in a variety of ways. First, planet-star flux ratios are a strong function of wavelength. 


owing to the presence of gaseous and aerosol absorption/extinction features (e.g., Marley 


et al. 1999, Sudarsky et ah 2003; Burrows et ah 2004, Cahoy et ah 2010). Thus, while a 


coronagraph may provide a raw contrast, or achieve a contrast floor, that can detect certain 
planet types at continuum wavelengths, this performance may not permit adequate SNRs in 
molecular absorption bands, thereby impeding characterization attempts. 


Second, as coronagraphs have inner and outer working angles that depend on wave¬ 
length, the variety of planets accessible for observation will depend on distance from the 























15 


Solar System. The interplay between planet-star flux ratio, wavelength, host star effective 
temperature, and inner and outer working angle constraints are depicted for different planet 
types in Figures |2]-|^ These contour diagrams are extremely useful for understanding the 
types of observations that could be made by a given mission architecture and, to our knowl¬ 
edge, are the first of their kind. Additionally, these diagrams depict an important reality 
of exoplanet spectral characterization; wavelength-dependent instrument inner and/or outer 
working angles can cause portions of a exoplanetary spectrum to be unobservable, depending 
on the apparent planet-star separation at the time of observation. 


Recall that we have defined a planet “type” based, in part, on insolation and size. Thus, 
Figures |2]-[^ show that a. 2X/D IWA for a 2-meter class telescope will strongly limit our 
ability to study warm Jupiters and Neptunes (at distances equivalent to 0.8 AU from the 
Sun) and Venus-like planets. Cool Jupiters and Neptunes (at distances equivalent to 2 AU 
from the Sun) are much more accessible. Interestingly, the presence of thick water clouds 
in the atmospheres of the cool Jupiters and Neptunes make these worlds very reflective 
compared to their warm counterparts (Sudarsky et ah 2000, Cahoy et al. 2010). Hence, the 
planet-star flux ratios in Figures and are broadly similar for these planet types, even 
though the insolation levels of the warm planet types is over six times larger. 


For our baseline parameters, integration times for cool Jupiters are relatively short, 
even out to 10 pc (Figure]^. Detecting the bottom of the 0.73 /rm methane band, which 
could indicate different planet formation scenarios as it can distinguish between different 
amounts of heavy element enhancement (Cahoy et al. 2010), can be accomplished with a 
10 hr integration time (or less) for distances as large as 7 pc and with 100 hr integration 
times out to 13 pc (Figure]^. Integration times required to detect the 0.73 /rm methane 
band are similar, and could be decreased by about 50% by observing both the 0.73 /im and 
0.89 /xm methane bands. Detecting the bottom of the 0.89 /im band in less than 100 hr is 
limited to cool Jupiters at distances of 6 pc or less, where the overall longer integration times 
for detecting the bottom of this band are driven up by falling detector quantum efficiencies 
at these wavelengths. Finally, these integration times are largely insensitive to raw contrast 
performance and exozodi levels (Figure 13). Regarding the latter, dust clearing by giant 


planets will likely prevent large exozodi levels in the vicinity of a giant (e.g., Papaloizou 


et al. 2007), limiting the range of realistic values for in Figure 13 


Resolution 70 spectra at constant SNR for Jupiter twins are limited by long integration 
times in methane bands at long wavelengths, which can exceed 10^ hr (Figures]^. Such long 
integration times are unrealistic (especially on a shared resource like WFIRST-AFTA—the 
Hubble Deep Field used less than 200 hr of exposure time). Only for continuum wavelengths 
between about 0.5-0.8 fim do integration times for SNR = 5 spectra fall below 10^ hr for a 
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Jupiter twin at 10 pc. While integration times for Jupiter twins around closer (i.e., less than 
about 5 pc) stars would be much shorter, telescope OWA constraints would limit study of 
twins in these nearby systems to cool stars (i.e., roughly 4,300 K or less for the orbit to be 
entirely within the OWA). Note, however, that the OWA constraint is less strict than the 
IWA constraint, as some fraction of the orbit, both towards gibbous and crescent phases. 


can occur inside the OWA, depending on orbit parameters and orientation (see, e.g., Greco 


& Burrows 2015). 


For the cool Neptunes, a coronagraph contrast noise floor of 10“^° would strongly limit 
observations at longer wavelengths (Figures]^, although detecting the 0.62 /rm methane band 
with integration times below 100 hr can be done out to distances of about 5 pc (Figure [l0|. 
As this band is relatively deep in the reflectance spectrum of cool Neptunes, integration 
times for detecting the bottom of the 0.62 pm band are longer, remaining below 100 hr only 
out to 3-4 pc. The 0.54 pm methane band is shallower, so detecting its bottom can be done 
with integration times comparable to the band detection times for the 0.62 pm band. 

Inner working angle constraints strongly limit the ability of a 2-meter class telescope 
to study Earths and super-Earths at 1 AU flux equivalent distances from their host stars. 
Covering the entire 0.4-1 pm spectral range at a distance of 5 pc (or more) requires solar 
twins (or hotter). For Earth twins, this implies planet-star flux ratios of 10“^°-10“^^, which 
could be improved by a factor of (i^p/i?®)^ for super-Earths (assuming similar reflectivity). 
Additionally, maintaining integration times below even 10^ hr at continuum wavelengths for 
these worlds at 5 pc would require degrading spectra to TZ = 10-20 to achieve SNR = 5 
(Figure]^, notably ruling out detection of the 0.76 pm oxygen A-band (Brandt & Spiegel 
|2014[ ). Such long integration times may be an unrealistic price to pay for characterizing 
potentially habitable worlds. Additionally, a,2\/D coronagraph IWA could interfere with 
our ability to detect the 0.94 pm water vapor band, a key indicator of habitability, for 5 pc 
Earths and super-Earths around Sun-like stars and cooler. These IWA constraints could be 


alleviated by investigating worlds at the outer edge of the habitable zone (Kasting et ah 


1993, Kopparapu et ah 2013[), but the low incident flux on these worlds pushes integration 


to times >10^ hr (assuming an Earth-like size and reflectivity), even at TZ = 10. 

Shorter integration times and more favorable contrast ratios argue for super-Earths 
as being more attractive targets than Earth-sized worlds. Even for these worlds, though, 
detecting the oxygen A-band with integration times less than 10^ hr requires targets to 


be nearer than about 3 pc (Figure 13). Only a few stars meet the distance and effective 
temperature requirements {a Centauri A and B, e Eridani), and most of these stars are in 
multi-star systems. Detecting the 0.94 pm water vapor band in super-Earth spectra with 
integration times below 10^ hr will be extremely challenging, if not impossible, owing in large 
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part to extremely low detector quantum efficiencies at these wavelengths. This outlook could 
be improved by using detectors with better quantum efficiency at red visible wavelengths, 
although these detectors tend to have larger dark currents, which might diminish any possible 
improvements. 

While our results show that characterizing Earth twins with a 2-meter class, coronagraph- 
equipped telescope will be unlikely, 10-meter class LUVOIR telescopes show much more 
promising results (Figure 16). Such a mission could achieve 77. = 70 visible wavelength 


spectra at SNR = 5 for Earth twins at a characteristic distance of 10 pc for integration 
times of order 100 hr. Falling stellar spectra and strong water vapor absorption drive up 
integration times at near-infrared wavelengths. Additionally, falling CCD quantum efficiency 
will drive up required integration times for species and band bottom detection at red visible 
wavelengths, especially for the 0.94 /im water vapor band. Using the near-infrared HgCdTe 
detector at these wavelengths can improve required integration times by an order of mag¬ 
nitude, even when accounting for the larger dark current and read noise in these detectors. 
Finally, note that telescope IWA constraints can prevent observations of large parts of the 
near-infrared wavelength range for Earth twins orbiting Sun-like stars at 10 pc, and exozodi 
levels larger than Nez = 10 begin to rapidly drive up integration times. 

Moving beyond requisite integration times and sensitivities, several of our assumptions 
about our instrumental, astrophysical, and planetary parameters warrant further discussion. 
Most notably, current experimental coronagraphs designed to achieve raw contrasts suitable 
for exoplanet characterization typically operate only over a ~20% bandpass (e.g., Trauger 


et al. 2013), which is limited by chromatic effects within the wavefront control system (al¬ 


though progress is being made on some aspects of this problem, e.g., Newman et al. 2015). 


Thus, if only a single coronagraph is used, observing over the full 0.4-1.0 /rm range used in 
this work would require 4-5 separate integrations, thereby driving up overall characterization 
time. This feature of coronagraphs is not entirely negative, as integration times in different 
bandpasses could be tailored to prevent achieving unnecessarily large SNRs in bright spectral 
regions. Also, while we include stellar leakage from the coronagraph, we did not consider how 
systematic errors from speckle subtraction might impact performance and increase required 
integration times, all of which would depend on telescope stability, the capabilities of the 


coronagraph wave-front control system, and how well speckles can be modeled (Krist et al. 


2008). 


Regarding astrophysical assumptions, recall that our noise models assumed that the star- 
planet systems were isolated in the instrument field-of-view. In reality, stellar companions 


(which are common, see Duchene & Kraus 2013) and background objects will contribute 


light to the observations, possibly within the planetary PSF. Stray light from a companion 
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could impact the systematic contrast noise floor and/or the raw contrast performance. While 
simplistic for companion star scenarios, our sensitivity studies to raw contrast performance 
{C) provide an estimate for how integration times would change if a companion were to 
increase the contrast background. Also, we assume smoothly varying exozodi structure that 
only depends on distance from the host star. Real exozodiacal disks will have more complex 
structure, such as clumps and bands, and studies of such structure would be a major science 
driver of a space-based coronagraph mission. 


While we have presented a range of planet models and types, the true diversity of ex¬ 
oplanets will be much greater. Additionally, some of our planet types (the warm and cool 
Jupiters and Neptunes) are based on models that do not include atmospheric photochem¬ 
istry. For the Solar System gas and ice giants, this chemistry leads to haze formation, which 
causes lower albedos at blue wavelengths (Karkoschka 1994). Thus, we are likely underesti¬ 
mating integration times for these worlds at the shortest wavelengths. Also, our placement of 
planet types at flux equivalent distances attempts, to first order, to maintain self-consistent 
atmospheric and cloud structures for these worlds. This treatment does not account for the 
shift of cooler stellar spectral energy distributions to longer wavelengths, where planets tend 
to be more absorbing ( Marley et al.||1999 ). So, our planet types around cooler stars would be 
absorbing more stellar flux, implying that an “absorbed flux” equivalent distance is further 
from the host star than our incident bolometric flux equivalent distance. 


5. Conclusions 

Space-based 2-meter class exoplanet characterization missions have the potential to 
provide spectra at moderate to high SNRs (i.e., 5 or better) for a variety of planets, some of 
which have no analog in our Solar System. We And that: 

• for cool Jupiters (i.e., Jupiters at 2 AU flux equivalent distance from the Sun), methane 
can be detected with integration times shorter than 100 hr out to a distance of 13 pc, 
and with integration times less than 10 hr at 5 pc. 

• spectra of cool Neptunes at 10 pc will be challenging (i.e., integration times approaching 
or exceeding 10^ hr), although methane could be detected in the spectra of these worlds 
with integration times below 100 hr out to a distance of 5 pc. 

• Earth twins are unlikely to be characterized, owing to long integration times and planet- 
star flux ratios smaller than 10“^° for a wide range of stellar effective temperatures. 

• super-Earths could be studied with low-resolution spectra and integration times below 
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10^ hr to a distance of 5 pc, but detecting features in higher-resolution spectra will be 
limited to worlds around a small handful of stars within 3 pc of the Sun. 

• falling CCD quantum efficiency at red visible wavelengths strongly limits observations 
of the 0.89 /xm methane band and the 0.94 /rm water vapor band. 


We note that the majority of the modeled observations of planet types presented here 
have Afexp oc Cn/Cp (to achieve a given SNR). Thus, major gains in driving down requisite 
integration times could be made by striving for larger instrument throughputs, better de¬ 
tector quantum efficiencies, and by using devices with low dark current (although achieving 
lower dark currents in modern detectors can lead to trade-offs with other detector proper¬ 
ties, Harding et ah 2015). Observational investigations will not have control over planetary 
properties, but we do have control over instrument performance requirements. 


A. Appendix: Noise Model Expressions 

This Appendix contains a brief discussion of our models of the planet count rate and 
count rates from key noise sources. 


A.l. Planet Count Rate 


The planetary signal is the incident stellar host spectrum, weighted by a wavelength- and 
phase-dependent planetary reflectance, and passed through the optics of a distant telescope. 
We assume that the host star has radius Rs and a spectrum of a blackbody, Bx (Tefj), where 
Teff is the stellar effective temperature. Thus, the specific flux density at a distance d from 
the stellar host is, 

F,^x{d) = TiBxins) 

Note that non-blackbody stellar spectra, either from observations or models, can be used 
for Fs^x{d )—we simply adopt blackbody spectra for efficiency, allowing us to rapidly model 
large swaths of parameter space. If the planetary target has radius Rp, orbits its host 
at a distance r, is observed at phase angle (i.e., star-planet-observer angle) a, and has a 
wavelength-dependent geometric albedo A, then the specihc flux density at a distance d 
from the planet is. 


Rr 


d 


-f = 7rA^ia)BxiRs) — -f 


Rs 


Rr 


d 


Fp,x{d) = A^a)F,,x{r) 


2 


(A2) 
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where <h(a) is the planetary phase function, which is generally a function of wavelength. 
Thus, the wavelength-dependent planet-star flux ratio is, 

which, of course, is only a function of the product of the geometric albedo and phase function 
(i.e., the phase-dependent planetary reflectivity), and the ratio of the planetary radius to 
the orbital separation. 

At the telescope, a distance d from the star-planet system, the planetary flux is passed 
through the optical components, with total throughput T (which can, in general, be wavelength- 
dependent) and spectral bandwidth AA, and converted into counts on the detector with 
quantum efficiency q (in dimensionless units of counts per photon). Thus, the planet count 
rate is given by, 

Cp = 7rg/par^Fp,A(d)AA , (A4) 

which has dimensions of counts per unit time. Here, h is the Planck constant, c is the speed 
of light, and /pa is a factor of order unity (computed using X, the photometric aperture 
length) that describes the fraction of light from the planet that falls within the photometric 
aperture. Inserting Equation |A2| and rearranging, we have 

= . (A5) 

Furthermore, for a spectrometer with constant resolution, TZ = A/AX, we can express the 
planet count rate as 

= . (A6) 


A.2. Zodiacal and Exozodiacal Light 


Local zodiacal light is known to be a function of ecliptic latitude and longitude (Levasseur 


Regourd & Dumont 1980), although Stark et ah (2014) demonstrated that a constant V- 


band surface brightness of = 23 mag arcsec ^ was a reasonably accurate representation. 
Thus, we take the wavelength-dependent zodiacal light count rate to be given by 


A (D 


i"0,A(l AU) 
Fo,v{l AU) 




(AT) 
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where Fq^x is the wavelength-dependent specific solar fiux density, -f©,!/ is the solar fiux 
density in V-band, hi is the photometry aperture area (with dimensions of arcsec^), and Fq v 
is the standard zero-magnitude V-band specific fiux density {Foy = 3.6 x 10“® W m“^ 

We use a square aperture with width XX/D along an edge, such that Q = {XX/D)'^, to most 
accurately represent a square grid of pixels (or spaxels) that would be used to isolate the 


planet. Note that the ratio Fe, x/FQy in Equation |A7| forces the color of the zodiacal light 
spectrum to match that of the Sun, which is a reasonable approximation, especially at large 
elongation angles from the Sun (Leinert et ah 1981). Inserting hi, and assuming constant 
spectral resolution, we have 


c, = TrgTX^ , 

^ AhcR FQy{l X\J) 


(AS) 


which is a strong function of wavelength. 


Following Stark et ah (2014), we define a “zodi” as the surface brightness of an exozodi- 


acal disk at 1 AU from a solar twin, which we label M^^y- Then the wavelength-dependent 
exozodiacal light count rate is given by. 


A {D 

c.. = .,T!JAA- - 


' X Auy F,^x{l AU) F,,v(l AU) 
r F,,v(l AU)F©y(l AU)' 




(A9) 


where is the number of exozodis in the disk. Note that the Fgy/FQy term ensures 
that exozodiacal light surface brightness scales with the intrinsic stellar brightness (at fixed 
orbital distance), and the (1 AU/r)^ term causes the disk surface brightness to decrease with 
increasing orbital separation according to the 1/r^ law. Again, inserting hi, and assuming 
constant spectral resolution, we have 


= irqTX^ 


X^ flAUV Fs,a(1 AU) F,y{l AU) 


Ve,Fo,ylO“^^^’^A.5 


(AlO) 


4hc7^ V r y Fs,v(l AU)F0,y(l AU)' 

For a discussion of the variety of different exozodiacal light treatments and exozodi definitions 


adopted throughout the literature, see [Roberge et ah (2012). Note that, while the zodiacal 


and exozodiacal light count rates are independent of telescope diameter, the ratio of the 
planetary signal to the zodiacal and exozodiacal light signals will depend on telescope size 
[Cp/ (Cz “1“ Cez) OC D ]. 


A.3. Leakage, Dark Counts, and Read Noise 


We consider several other key sources of noise. Stellar light leaked through the coron- 
agraph (“leakage”) can contribute a large number of noise counts as the host star is many 
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orders of magnitude brighter than the target planet. The leakage count rate is determined 
from the stellar photon flux passed through the observing system, diminished by C (the 
coronagraph design raw contrast), and is given by 

(f J = ^ (^J . (All) 

where the hnal step assumes constant spectral resolution. Detections are possible for sources 
fainter than the raw contrast, but coronagraphs will have systematic noise floors that limit 
how far below the raw contrast observations can go. The noise floor limit is due to residuals 
in the subtracted stellar PSF, and is typically a factor of ten below the raw contrast. This 
limit is important to keep in mind when discussing small, low-albedo, and/or large separation 
planetary companions. 

For a given dark current Dg-, in counts per pixel per unit time, the dark count rate is 


cd Dg— fVpjx , 


(A12) 


where iVpix is the number of detector pixels contributing to the spectral element. For imaging, 
the contributing pixels are those that fall within the aperture (D), so that, if a pixel has an 
angular diameter dpix, then 


N ■ ■ = 

’ plX,l 


4D 


(A13) 


where the sub-script ‘i’ has been introduced for ‘imaging’. If the pixel diameter is determined 
by the diffraction limit at some wavelength Aq, giving 6'pix — Xo/2D, then the number of 
contributing pixels is simply 


AT. . = 1^ 

’plX.l 

TT 



(AM) 


Determining the number of contributing pixels for a spectrometer is more complicated. 
For an integral held spectrometer (IFS), angularly-resolved spectroscopic observations are 
obtained by placing a grid of small lenses (called “lenslets”) in the focal plane, spectrally 
dispersing the light from each lenslet, and recording the resulting spectra at the detector. 
For this setup, the number of contributing pixels is 

4D 

ApixjS ^pixAA —A , (Al5) 

■^^lens 

where ‘s’ is for ‘spectroscopy’, Upix is the number of pixels per unit wavelength designed 
for each lenslet spectrum, and ^lens is the angular diameter of an individual lenslet (so that 
is the number of lenslets over which the aperture is distributed). For a spectrum 
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sampled at two pixels per spectral resolution element, and spread over AA^hpix pixels in the 
horizontal/spatial dimension, we have 


27^ 


‘^pix 


''hpix 


2AA^hpix 

AX 


(A16) 


which gives 

Api,,, = SAAhpi,^^ . (A17) 

^^lens 

Again, if the lenslets are sized to yield the diffraction limit at some reference wavelength, 
then we would have 

32 / XA\ ^ 

Api,,, = AAhpi.- — . (A18) 

TT \ Ao / 

Note that the lenslet array typically only covers the inner portion of the overall imaging 
area, with coverage limited by detector space for recording individual lenslet spectra. Thus, 
lenslet arrays and, more importantly, their associated detectors can set a tighter constraint 
on the OWA than the coronagraph. 


Read noise is computed using the number of reads per observation, Aread, the number 
of contributing pixels, and the read noise counts per pixel. Re-, and is given by 


cr 


A^pixA^read 


Re- , 


exp 


(A19) 


where dividing by the exposure time, Atexp, permits expressing read noise counts as a rate. 
The number of reads is determined by assuming a maximum exposure time, Atmax (often 
limited by cosmic ray strikes), with 


Atexp 

’ rpptn 


AA 


so that Equation A19 can be written as 

Cr = 


A ■ 

pix p 

JLp 


At 


(A20) 


(A21) 


Note that integer quantities (e.g., numbers of pixels) can be determined by rounding the 
expressions above, which are, for simplicity, expressed as non-integer quantities. 

For completeness, internal thermal noise can be included and has a count rate given by 
Cth = 7TqQ^eeysBx{Tsys) AX = nqesysBxiRys) ’ (^22) 
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where egys is the effective emissivity for the observing system (of order unity), and T^yg is the 
system temperature. As we focus primarily on the visible wavelength range throughout this 
paper, internal thermal contributions are negligible, so we ignore them throughout. 
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B. Tables and Figures 


This preprint was prepared with the AAS IATeX macros v5.2. 
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Table 1. Symbol Usage 


Symbol Description 

A wavelength-dependent planetary geometric albedo 

a. planet phase angle 

B\ Planck function 

C coronagraph design contrast 

C’noise total number of noise counts on detector for a spectral element 
Cb total number of background counts on detector for a spectral element 

Ctot total number of counts on detector for a spectral element 

Cb total background count rate on detector for a spectral element 

Ctot total count rate on detector for a spectral element 

cd dark current count rate on detector 
Cez exozodiacal light count rate on detector 
Cp planetary count rate on detector 

cr read noise count rate on detector 

cik leakage count rate on detector 

Cz zodiacal light count rate on detector 

c speed of light 

D telescope diameter 

Dg- dark current 

d distance to observed star-planet system 

Fp^x planetary specific flux 

Fq X stellar specific flux 

^0.A solar specific flux 

Fe,v V-band solar specific flux 

Fs,v V-band stellar specific flux 

Fq y standard zero-magnitude V-band specific flux 

/pa fraction of planetary light that falls within photometric aperture 
h Planck constant 

A wavelength 

AA spectral element width 

V V-band zodiacal light surface brightness 

^ez,v V-band exozodiacal light surface brightness 

Nez number of exozodis in exoplanetary disk 

F^read number of detector reads per observation 
-^pix,i/s number of contributing pixels for imaging/spectroscopy 
^^hpix number of horizontal/spatial pixels for dispersed spectrum 
^pix pixels per unit wavelength at detector for each lenslet spectrum 
Q photometry aperture size, expressed as A{XX/D)^ 

^>(q:) wavelength- and phase-dependent planetary phase function 

q detector quantum efficiency 

7Z instrument spectral resolution 

Rp planetary radius 

Rs stellar radius 

R^- read noise counts per pixel 

r planet-star distance 

SNR signal-to-noise ratio 

T telescope and instrument throughput 

^IWA coronagraph inner working angle 

^OWA coronagraph outer working angle 

0pix detector pixel angular diameter 

^lens lenslet angular diameter 

Teff stellar effective temperature 

Atexp exposure time 

Atrnax detector maximum exposure time 

X width of photometric aperture, as multiple of A/D 
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Table 2. Baseline Astrophysical Parameter Values 


Parameter 

Description 

Adopted Valne 

a 

planet phase angle 

90° 

Foy 

standard zero-magnitnde V-band specific flnx 

3.63 X 10“® W 

Mzy 

V-band zodiacal light snrface brightness 

23 mag arcsec”^ 

Mezy 

V-band exozodiacal light snrface brightness 

22 mag arcsec“^ 

Ne. 

nnmber of exozodis in exoplanetary disk 

1 


Table 3. Baseline Telescope and Instrument Parameter Values 


Parameter 

Description 

Adopted Valne 

C 

coronagraph design contrast 

10-9 

De- 

dark cnrrent 

5 X 10“^ sec-i 

D 

telescope diameter 

2 m 

/pa 

fraction of planetary light that falls within photometric apertnre 

0.87 

A/Vhpix 

nnmber of horizontal/spatial pixels for dispersed spectrnm 

3 

q 

detector qnantnm efficiency 

Eqnation 

7^ 

instrnment spectral resolntion 

70 

Pc- 

read noise connts per pixel 

0.1 

T 

telescope and instrnment thronghpnt 

0.05 

SiwA (A/D) 

coronagraph inner working angle 

2 

^'owA (A/D) 

coronagraph onter working angle 

10 

^^max 

detector maximnm exposnre time 

1 hr 

A 

width of photometric apertnre, as mnltiple of A/D 

1.5 
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Fig. 1.— Scaled opacities for different species and aerosols between 0.4-1.0 /ini aX IZ = 70, 
demonstrating the nnmerous bands and featnres present across the visible wavelength range. 
Scalings for gaseous absorption are constant with wavelength for each species and arbitrary; 
these were chosen to highlight relevant bands. The haze opacity is for extinction by a Titan 
tholin-lke haze, and the Rayleigh scattering opacities are scaled to that of CO 2 at 0.4 pm 
(and divided by 10). 
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Fig. 2.— Planet-star flux ratio color contours for a variety of Jupiter-sized worlds at fixed flux 
equivalent distances around main sequence stars of different effective temperatures. Planet 
types are: (a) a gray planet with A^{a) = O.S/tt and Ftoa = 50.5 W m“^ (top-left), (b) a 


Jupiter twin with Ftoa = 50.5 W m ^ and reflectivity from Karkoschka 

(1 

998 

(top-right). 

(c) a cool Jupiter with Ftoa = 342 W m ^ and reflectivity from 

Cahoy et a . 

(2010 

) (bottom- 

left), and (d) a warm Jupiter with Ftoa = 2140 W m ^ and reflectivity from 

Cahoy et ah 


(2010) (bottom-right). Inner and outer working angle constraints (at quadrature; dashed), 
at 2X/D and lOA/H, are shown for a 2-meter class telescope and assuming the planet-star 
system are at 5 pc (yellow) and 15 pc (orange). 
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Fig. 3.— Same as Figure but for Neptune-sized worlds. Planet types are: (a) a Neptune 


twin with Ftoa = 1.5 W m ^ and reflectivity from 

Karkoschka 

(1998) (to 

>left), (b) a cool 
) (top-right), and 

Neptune with Ftoa = 342 W m ^ and reflectivity from 

Cahoy et al. (2010 

(c) a warm Neptune with Ftoa = 2140 W m ^ and reflectivity from ICahoy et al. (2010) 


(bottom-left). Inner and outer working angle constraints are also as in Figure]^ 
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Fig. 4.— Same as Figure but for terrestrial worlds. Planet types are: (a) a Venus twin 
with Ftoa = 2610 W m~^ and reflectivity from an application of a line-by-line radiative 
transfer model ( Meadows &: Cri^|1996 ) to Venus ( Arney et al.||2014 ) (top-left) (b) an Earth 
twin with Ftoa = 1370 W m“^ and reflectivity from Robinson et ah (2010) (top-right), 


and (c) a super-Earth with l.SR®, Ftoa = 1370 W m ^ and Earth-like reflectivity from 


Robinson et ah (2010) (bottom-left). Inner and outer working angle constraints are also as 


in Figure 
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Fig. 5.— Contours of stellar effective temperature around which a certain planet type 
achieves a given planet-star flux ratio. Jupiter-sized cases are the same as Figure while 
the Earth case is as in Figure Inner and outer working angle constraints are also shown, 
and follow those in Figure 
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Fig. 6.— Contours of integration time required to achieve SNR = 5 for a Jupiter twin (left) 
and a cool Jupiter (right) at 10 pc. These planet types are placed at the flux equivalent 
distance of main sequence stars of different effective temperatures, as in Figure Note the 
two different color-bar scales. Regions excluded/limited by a 2X/D and lOX/D inner and 
outer working angle (respectively) are demarcated by a dashed line. 
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Fig. 7.— Same as Figure but for a cool Neptune placed at the flux equivalent distance, 
as in Figure]^ Regions with planet-star flux ratios smaller than 10“^° are above the labeled 
dotted line. 
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Fig. 8.— Same as Figure but for an Earth twin (left) and a 1.5i?0 super-Earth (right) 
at 5 pc. Both are placed at Earth’s flux equivalent distance (as in Figure]^. To maintain 
reasonable integration times, higher-resolution spectra {TZ = 70) have been degraded to 
77 = 10 and 77 = 20 for the Earth twins and super-Earths, respectively. Regions with 
planet-star flux ratios smaller than 10“^° are above the labeled dotted line. 
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Fig. 9.— Contours of integration time required to achieve SNR = 5 in the bottom of 
key absorption bands as well as contours of required integration time to detect a band at 
SNRband = 5 for cool Jupiters. We highlight the 0.73 /im and 0.89 /xm methane bands. The 
lower right plot shows how the band detection time decreases when both methane features 
are used. Regions excluded/limited by a 2\/D and lOA/D inner and outer working angle 
(respectively) are demarcated by a dashed line. 
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Fig. 10.— Similar to Figure except now for cool Neptunes. Here, for detecting the bottom 
of features, we highlight the 0.54 pm and 0.62 pm methane bands, as the longer wavelength 
bands have planet-star flux ratios of 10“^° or smaller in their bases. For band detections, 
we use the 0.73 pm methane band, as its depth leads to shorter required integration times 
than bands at shorter wavelengths. The lower right plot shows how the band detection time 
decreases when both the 0.62 pm and 0.73 pm methane features are used. 
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Fig. 11.— Similar to Figure]^ except now for Earth twins. Here, we focus on the 0.76 pm 
A-band of molecular oxygen and the 0.94 pm water vapor band. 
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Fig. 12.— The same as Figure [TT| except for 1.5i?0 super-Earths. 
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Fig. 13.— Contours of integration time required to achieve SNR = 5 in the bottom of the 
0.73 /im methane band (top) or to detect this band at SNRband = 5 (bottom) for a cool 
Jupiter at a distance of 10 pc for various levels of exozodis and raw contrast. Note that the 
IWA limits the range of that can be investigated (as the wavelength and distance are 
hxed). 
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Fig. 14.— Similar to Figure [T^ except for cool Neptunes at a distance of 5 pc. Here the 
0.62 pm methane band is used for the feature bottom values, while the 0.73 pm methane 
band is used for the band detection values. 
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Fig. 15.— Similar to Figure [T^ except for the required band detection time for the 0.76 /xm 
A-band of molecular oxygen for super-Earths at 3.7 pc. 
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Fig. 16.— Results for a 10-meter class space-based telescope and coronagraph, capable of 
10“^° raw contrast between an IWA of 3X/D and an OWA of 20X/D. For Earth twins 
around stars of different effective temperatures, panels show the integration times to achieve 
SNR = 5 and sensitivity to exozodi levels for integration times in the A-band of molecular 
oxygen at 0.76 /xm. 
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Fig. 17.— Simulated spectral observation of an Earth twin around a solar twin at 10 pc for a 
LUVOIR mission assuming a 200 hr integration. The dark-red line is a noise free spectrum, 
while the yellow points are a simulated spectrum (with l-cr error bars in gray). 











































































